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Several recent helicopter vibration reduction 
research programs of the United Technologies 
Research Center (UTRC) are described. Results of 
analytical investigations to develop, understand, 
and evaluate potential helicopter vibration reduc- 
tion concepts are presented in the following areas: 
identification of the fundamental sources of vibra- 
tory loads, blade design for low vibration, appli- 
cation of design optimization techniques, active 
higher harmonic control, blade appended aeromechan- 
ics! devices, and the prediction of vibratory air- 
loads. Primary sources of vibration are identified 
for a selected four-bladed articulated rotor oper- 
ating in high speed level flight. The application 
of analytical design procedures and optimization 
tecnniques are shown to have the potential for 
establishing reduced vibration blade designs 
through variations in blade mass and stiffness 
distributions, and chordwise center-of-$ravity 
location. Analytical evaluation of a computerized 
generic active controller for implementing higher 
harmonic control indicates the potential for good 
controller performance and extensive fuselage vi- 
bration reduction with low pitch amplitudes for 
three controller approaches investigated. Explora- 
tory evaluation of a passive tuned blade tab con- 
cept indicates considerable sensitivity of vibra- 
tory load alleviation to design parameters with an 
improvement in inplane hub excitation but an 
increase in vertical excitation. The prediction of 
vibratory airloads, attributable to rotor/wake, 
rotor/fuselage, and rotor/erapennage interactional 
aerodynamics, is al^o described. 

Introduction 

For future helicopters to reach their full 
potential, significant reduction in vibration must 
be accomplished. Helicopter vibration is becoming 
an increasingly important consideration because of 
requirements for crew and passenger comfort as well 
as increased reliability of structural components 
and on-board electronic equipment. A vibration- 
free weapons platform is also a requirement for 
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some military helicopters. As helicopters are 
required to fly more at both faster cruise speeds 
and slow "t rans it ion M speeds for nap-of-the-earth 
flying, the need to minimize vibration becomes more 
important, and it becomes necessary to consider 
alternative and complementary approaches for vibra- 
tion reduction. 

Helicopter vibration research has been a 
primary activit} at the United Technologies 
Research Center (UTRC) £ past six years. 

Various helicopter vibration related programs 
involving analytical investigation have been con- 
ducted in the following areas: 

• Blade design for low vibration 

• Design optimization techniques applicable 
to vibration reduction 

• Active higher harmonic control for vibra- 
tion alleviation 

• Blade appended davices >r vibration alle- 


• Prediction of vibratory airloads (rotor, 
fuselage, empennage). 

In the Ref. I analytical investigation, rotor 
vibratory response and loads transmitted to the 
fuselage were predicted and analyzed to determine 
the relative contributions and sources of the 
various components of blade force excitation. 
Primary sources of vibration were identified for a 
selected four-bladed articulated rotor operating in 
high speed level flight. Subsequently, blade modal 
shaping (Ref. 2), frequency placement, structural 
and aerodynamic coupling, and intermodal cancella- 
tion were investigated to systematically identify 
and evaluate blade design parameters that influence 
vibratory airloads, blade modal response, hub 
loads, and fuselage vibration. Through variations 
in blade mass distribution, stiffness distribution 
and chordwise center-of-gravity location, blade 
designs were developed with predicted reductions in 
vibration. These designs remain to be validated by 
test. 

An automated optimization procedure is being 
developed at UTRC for the rotor blade design 
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process. A simplified approach for minimizing 
vibration (Refs. 1, 2) has been developed and 
applied ^Ref. 3). A modal analyris was used to 
calculate key vibration parameters. Through 
frequency placement and modal shaping with a 
constrained optimization program, COPES/CONMIN, 
(Refs. 4, 5) a blade design was determined for 
reduced fuselage vibration. A forced response 
aeroelaetic analysis (Ref. 6) was used in the 
process to identify cv modes, the desired 
frequency placement and modal shaping criteria, and 
to perform a final calculation of the vibration 
characteristics of the new blade designs. 

A computerized generic active controller was 
developed for alleviating helicopter vibration by 
closed-loop implementation of higher harmonic 
control (HCC) (Refs. 7, 8). This controller 
provides th<* unique capability to readily define 
and evaluate many different algorithms by selecting 
from three controller approaches (deterministic, 
cautious, and dual), two linear system models 
(local and global), and several methods of limiting 
control. A non-linear aeroelastic rotor analysis 
(Ref. 6) was used to evaluate alternative control- 
ler configurations as applied to a fo . - bladed H-34 
rotor mounted on the NASA-Amea Rotor Test Apparatus 
used to represent the fuselage. It will be shown 
that excellent controller performance was predicted 
for all three controller approaches for steady 
flight conditions having moderate to high values of 
forward velocity and rotor thrust. Reductions in 
vibration from 75 to 95 percent were predicted with 
HHC pitch amplitudes of less than .ie degree. Good 
transient vibratio alleviation '.as also predi ed 
for short duration maneuvers involving a sudden 
change in collective pitch. 

Analytical evaluations of aeroelastic devices 
appended lo helicopter rotor blades have been 
conducted to determine their pocen.ial for reducing 
hub shears and vibratory control loads (Ref. 9). 
The results for a passive tuned tab shall be 
discussed. 

The developm nt of methodology to pt edict 
vibratory blade i irioads has proceeded at UTRC 
along with the determination of primary airload 
sources. Wake, airflow, and airload methodology 
have advanced for both low and high speed flight 
(Refs. 10-19). For example, a fint level general* 
ization of the forward flight rotor wake has 
recently been formulated for use in unsteady air- 
load calculations (Refs. 10, 11, 12). Synthesiza- 
tion procedures for incorporating unsteady airfoil 
teat data in rotor aeroelastic response methods 
have been developed (Ref. 18). Also, methodology 
for predicting vibratory airload excitation at tail 
surfaces due to the rotor; wake has been developed 
and initial validation has been performed (Ref. 
19). The aerodynamic interaction of the fuselage 
on the rotor vibratory airloads has been analytic- 
ally demonstrated (Ref. 17). Development of a 
computer method for predicting the induced unsteady 
vibratory excitation of ttu rotor wake on the fuse- 


lage ia in progress. 

The aforementioned investigations are 
described in the following sections. 

Blade Design for Vibration Reduction 

An analytical investigation was conducted to 
develop an understanding of the importance and role 
played by blade design parameters i\i rotor vibra- 
tory response and to design an advanced blade for 
reduc vibration based upon thiL understanding. 
This investigation was conducted at UTRC by Taylov 
j.« M, Various design approaches were examined 
for a fc --bladed articulated rotor operating at a 
high-speed, level flight condition. Blade modal 
shaping, frequency placement, structural and aero- 
dynamic coupling, and intermodal cancellation were 
investigated to systematically identify and 
evaluate blade design parameters that influence 
blade vibratory airloads, blade modal responses, 
hub loads, and fuselage vibration. 

The baseline rotor system selected for the 
investigation wae a four-bladed articulated rotor 
system, similar in design to the Sikorsky S-76 
rotor system, but without vibration alleviation 
devices ami tip sweep. Tip sweep can provide vi- 
bration alleviation through aeroelastic coupling 
between the blade flatwise and blade torsion modes. 
However, It was decided to omit tip sweep from the 
baseline blade design to study the potential of 
improving the vibration characteristics of a basic 
rectangular planform blade. A complete set of 
baseline blade properties is presented in Ref. 1. 
The blade spanwise mass distribution for the base- 
line and modified blade designs will be presented 
in the following section on blade design optimiza- 
tion. Some of the relevant properties of the base- 
line blade are the blade twist (-10 deg, nonlir- 
ear), blade weight (100 lb) and non-dimensional ized 
natural frequencies (flatwise modest rigid jody ■ 
1.03 per rev, first elastic ■ ' .15 per rev, second 
elastic "4.9 per rev; edgewise modes: rigid body 
■ 0.26 per rev, first flexible " 4.7 per rev; 
torsion mode: first elastic * 5.30 per rev). The 
baseline blade pitch axis, elastic axis, aerodynam- 
ic center, ana center-of-gravity (CG) were nomin- 
ally located at the 25 percent chord (outboard, CG 
was at 26Z). The flight condition selected was a 
high speed cruise condition of 160 kt. A 10,800 lb 
rotor lift and an 1190 lb propulsive force were 

selected to be representative for a helicopter the 
size of the S-76. 

The enalyticel simulation used for this etudy 
was the G400 analysis documanted in Ref. 6. This 
computer analysis is based on the Celerkin method 
nd uses precalculated uncoupled normal blade 
modes. A time history solutior airloads, blade 
responses, end vibratory blade r.j’o forces and mo- 
ments due to the rotor ere calculated baaed on 
modal coupling within the analysis. For this 

study, the fuselage vibrations were determined by 
means of a measured 3-76 mobility transfar matrix 
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between the hub loads *nd vibrations at several 
locations in the airframe. This approach was used 
to separate rotor design effects from hub impedance 
and airframe dynamic effects in the analysis. 
Uniform inflow was used in this in.tial exploratory 
study to model the rotor downwarh for the high 
speed flight condition. The influence of variable 
inflow on vibratory airloads will be included in 
subsequent blade design studies. 

In order design blades for reduced heli- 
copter vibration, it is necessary to first acquire 
a fundamental understanding of the interrelation- 
ship of the fuselage vibrations, the hub loads in 
the fixed system, the blade root loads in the 
rotating system, the blade vibratory response, and 
the blad^i airloads. It is then necessary to 
identify the primary contributors of th^ various 
harmonic modal comf^oents to vibratory excitation. 
As indicated in /ig. 1, for a helicopter with a 
four-bladed articulated rotor, the predominant 4P 
(4 per rev) cockpit and cabin vibrations are mainly 
produced by the 4P vertical and inplane (lateral 
anc longitudinal) hub forces, since for a low-off- 
set articulated rotor the 4P moment contribution is 
of secondary importance. For a four-bladed rotor 
in a steady-state condition, the 4P hub fixed sys- 
tem vesical shear results directly from the azi- 
muthal summation of the 4P blade root rotating 
vertical shear from all blades. When the contri- 
butions of all four blades are sunmed in the fixed 
system, the 4P hnb fixed system inplane shears 
result from the summation of 3P and 5P blade root 
rotating inplane shears. As indicated in Fig. 1, 
the contribution of the 5P inplane shears were 
found early in this investigation to be smal 1 rela- 
tive to the 3P contribution. 



Fig. 1. Primary Contributors to Helicopter Vibra- 
tions (4-Bladed Articulated Rotor, High 
Speed Flight) 


With the awareness of the above information at 
the beginning of the investigation, the analytical 
study was directed toward identifying the primary 
blade response contributors and related airload 
harmonics, as will be discussed. With this infor- 
mation and an understanding acquired Irom sensitiv- 
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ity studies, several blades were designed and 
reported in Ref. 1 which resulted in lower pre- 
dicted vibration levels. Results for the baseline 
and the two final modified blade designs (Designs A 
and B) are summarized herein. It should be noted 
that these blade designs were achieved via closed- 
loop optimization techniques as discussed in the 
next section. 

The predicted 4P hub loads and fuselage vibra- 
tions for the baseline and two modified blade 
designs are presented in Figs. 2 and 3, respec- 
tively. For the baseline articulated rotor, the 
vibratory hub load components with the most 
influence on fuselage vibration are, in order of 
importance, the longitudinal, lateral, and verticil 
shears. On other helicopters having a different 
mobility matrix, the vertical shear may be more 
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Fig. 2. Predicted Vibratory Fixed-System Hub 
Loads for the Baseline and Modified Blade 
Designs (V * 160 kt) 
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Fig. 3, Predicted 4 Per Rev Fuselage Vibrations 
for t\e Baseline and Modified Blade 
Designs (V ■ 160 kt) 
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important than the inplane shears. For hingeless 
or be~ringless rotors, the hub moments can become 
important. Of the fuselage vibrations shown in 
Fig. 3, the dominant vibrations for the aircraft 
studied are the heel-slide vibrations for the pilot 
and copilot. Important vibrations are those af- 
fecting crew and passenger comfort. Thus, it was 
considered a good objective to reduce the cockpit 
seats and cabin vibrations to approximately 0.1 g 
while reducing the heel-slide vibrations to the 
extent possible. As shown in Fig. 3, the predicted 
fuselage vibrations for modiried Designs A and B 
are lower by about 50 percent (except the small 
pilot seat vertical vibration) compared to the 
baseline vibrations. Of particular importance are 
the seat and cabin vibrations which are approaching 
the 0.1 g level. These results indicate that ac- 
ceptable seat and cabin vibration levels can be 
attained, at least theoretically, by improved blade 
design and without the use of vibration treatment 
equipment (bifilars, absorbers, etc.). The vibra- 
tory hub loads that produced these vibrations are 
shown in Fig. 2 to also have been reduced signif- 
icantly. 

To achieve the vibration levels predicted 
for Designs A and B, the following modifications 
were made to the baseline blade. For Design A, a 
combination of blade spanwise mass redistribution 
and increased edgewise blade stiffness was used to 
change the blade mode shapes and increase the un- 
coupled frequencies of the first elastic flatwise 
and edgewise modes from 2.75 to 3.4 per rev and 4.8 
to 5.8 per rev, respectively. Modal shaping tech- 
niques described in Refs. 1 and 2 and summarised in 
the next section on optimization were applied in 
addition to frequency placement techniques, A 
small change in blade weight (100 to 104 lb) re- 
sulted. Increasing the flatwise frequency r.o 3.4 
oer rev improved the inplane shear, but did not 
significantly reduce the vertical shear, as shown 
in Fig. 2. For Design B, the generalized airload 
producing a large response of the first flatwise 
mode was decreased by moving forward the CG loca- 
tion of the outer 20 percent of the blade from 26 
tc 24 percent of the chord. The basis for these 
design changes are discussed below. 

In order to understand the source of the 
rotating blade root shears, the G400 analysis was 
used to decompose the vertical and inplane shear 
components (edgewise and radial) into the indi- 
vidual contributing components. This is exempli- 
fied in Fig. 4 for the 3P rotating edgewise shear 
transmitted to the hub by the baseline bladj. The 
force vectors are presented in polar format showing 
the amplitude and phase of the total shear and 
contributing components. It is important to note 
that the two principal components due to rigid body 
motion (rigid flapping coriolis and lag inertia) 
vectorially combine with the drag component to form 
a small contribution to the total. 

As a result of this cancellation effect, the 
elastic blade components of shear acquire a primary 
role. In particular, the elastic flatwise inertia 


component is the predominant contributor to the 3P 
edgewise shear. This component arises from the 
flatwise mode coupling with the edgewise mode 
associated with blade twist and collective pitch. 
In Fig. 5, it is shown that the elastic flatwise 
contribution is also important to the radial shear 
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Fig. 4. Predicted 3 Per Rev Rotating Blade Root 
Edgewise Shear for the Baseline Blade (V 
- 160 kt) 
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Fig. 5. Predicted 3 Per Rev Rotating Blade Root 
Radial Shear for the Baseline Blade (V ■ 
160 kt) 
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component. It thus became evident that significant 
intermodal coupling occurs, and reduction of the 
vibratory inplane hub loads is heavily dependent 
upon reduction in the response of the first elastic 
flatwise mode in addition to the first elastic 
edgewise mode as summarized in Fig. 1. It was also 
noted that significant vectorial cancellation of 
components of shear due to phase as well as ampli- 
tude differences occurred. Thus, the predicted 
total shears are very sensitive to the accuracy of 
the simulation analysis. 

Also, it is described in Ref. 1 that signi- 
ficant interharraonic and intermodal coupling occurs 
between the blade flatwise response and the air- 
loads. This is indicated in Fig. 1 as a "cascade 
effect" in which modal motions at one harmonic 
induce modal motions at other harmonics through the 
coupling with the harmonics of airloads. For for- 
ward flight, the one per rev airloads created by 
flapping (or cyclic pitch), required for aircraft 
trim, result In one per rev response of all blade 
modes. The resulting one per rev motions create 2, 
3, 4, and 5 per rev response of the blade modes and 
so on. The end result for the rotor investigated 
was that significant 3P and 4P airloads were gener- 
ated which excited the important 3P and 4P 
responses shown in Fig. 1. The rotor thus largely 
excites itself through the blade motion-airload 
cascade effect. Also, as will be described, har- 
monic inflow and airloads due to rotor wake effects 
provide vibratory excitation which, for steady 
level flight, is large at low speeds and less at 
high speeds. 

For modified Design A, the r* duct ion of the 3P 
root edgewise shear and its r . onents is shown in 
Fig. 6 (note the scale change relative to Fig. 4). 
The important elastic inertia components were sub- 
stantially reduced. The reduction in the predicted 
3P and 4P modal responses are shown in Fig. 7. 
Although the inplane hub shear was largely reduced 
for Design A, the vertical hub shear was not, as 
shown in Fig. 2. This was attributed to an over 
reduction of the second elastic flatwise mode con- 
tribution and the fact that the second mode pro- 
vided phase cancellation with the first mode. 
Although cancellation between the two modal contri- 
butions to vertical shear could have been pursued 
by blade design retuning with mass and stiffness 
variations, it was decided not to depend upon this 
cancellation due to intermodal phasing which could 
change with rotor configuration or flight condi- 
tion. Instead, the 4P vertical shear was reduced, 
as reported in Ref. 1, by changing the blade out- 
board center-of-gravity location to reduce the 
generalized airload. 

Changing the center-of-gravity (CG) location, 
over the outer blade region, from 26 to 24 percent 
of the chord moved the CG forward of the elastic 
axis (EA) and influenced the blade airloads through 
inertial couplirg between the blade flatwise and 
torsion responses. The predicted change in the 
torsion response time history around the azimuth 
for this modified blade Design B is shown in Fig. 
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Fig. 6. Predicted 3 Per Rev Rotating Blade Root 
Edgewise Shear for Modified Blade Design 
A (V • 160 kt) 
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Fig. 7. Predicted 3P and 4P Modal Response for 
the Baseline and Modified Blade Designs 
(V * 160 kt) 










8. 7jr the baseline design, the torsion response 
is characterised by a strong nose-down deflection 
on the rotor advancing side. For blade Design B, 
the nose-down response on the advancing side was 
eliminated leaving a waveform that consists of a IP 
response, maximum nose down at 180 degrees azimuth, 
plus smaller contributions of higher harmonics. 
The mechanism involved in this phenomenon is the 
CG-EA offset. When the blade tip bends downward on 
the advancing side, a forward CG induces a nose-up 
torsion increase on the blade near the blade tip. 
This moment counteracts the nose-down aerodynamic 
pitching moment due to high Mach number on the 
advancing side. It was found that, when the CG is 
moved further forward, the increased advancing side 
nose-up response can produce increased vibration. 
Blade tip sweep is an alternate approach that has 
been used to reduce the advancing side vibratory 
airloads through blade aeroelastic torsional defor- 
mation. 



Fig. 8. Predicted Effect of Outboard CG Location 
on Blade Torsion Mode Response (V * 160 
kt) 

The reduction of the torsion mode response of 
modified Design B i-* shown in Figs. 7 and 8. The 
reduction of the 4P vertical shear relative to the 
baseline and Design A blade levels is shown in Fig. 
2. Relative to Design A levels, a small further 
reduction in 4P fuselage vibration is indicated at 
all fuselage locations in Fig. 3. The resulting 
fuselage ’ibratioi are generally about half the 
baseline blade levels. 

In summary, the predicted results for the 
modified blade designs in Ref. 1 indicate the po- 
tential for su’ stantially reducing helicopter vi- 
bration by viable changes in blade design. Further 
study in the area of blade design for vibration 


could possibly result in rotors with vibration 
levels that are acceptable without the use of 
vibration alleviation devices or higher harmonic 
control. However, it is recognized that the pre- 
dicted results have not yet been substantiated with 
test. Due to the absence of a systematic, consis- 
tent set of experimental data for vibratory blade, 
hub, and fuselage loads, the G400 analysis, like 
other rotcr aeroelastic analyses, has not been 
validated for vibration prediction. The high 
degree of interharmonic and intermodal coupling and 
the predicted cancellation effects of large compon- 
ents of hub loads make the results sensitive to the 
prediction accuracy of the analysis for the indivi- 
dual components. In fact, differences ot the pre- 
dictions with the vibration results of preliminary 
exploratory tests at Sikorsky have been noted. The 
provision of a model test rig to systematically and 
accurately measure vibratory hub loads is currently 
being pursued at UTRC to validate computer codes 
and evaluate new blade designs for vibration. 
Also, application of blade design techniques for 
vibration to other rotor types (hingeless, bearing- 
less) is underway and application to other flight 
conditions (particularly low speed) with variable 
inflow is planned. The inclusion of design opti- 
mization techniques has been initiated as described 
below. 

Rotor Blade Design Optimization for Vibration 

The optimization approach discussed below is 
part of an ongoing effort at UTRC to develop a 
general automated procedure for rotor blade design. 
This procedure can be used to determine the 
necessary geometric, structural, and material 
properties of a rotor system to achieve desired 
design objectives relating to vibration, stress, 
and aerodynamic performance. This section concen- 
trates on the approach used for helicopter vibra- 
tion summarized in Ref. 3. Based on the analytical 
studies discussed above, a simplified vibration 
analysis has been developed for use in conjunction 
with a forced response analysis in the optimization 
process. TnU simplified analysis significantly 
improves the efficiency of the design process. 

Optimization Approach 

As shown in Fig. 9, the approach for rotor 
blade design has been formulated as three separate 
component optimization problems concerned with 
areas such as vibration, stress, and aerodynamic 
performance. Appropriate constraint functions are 
formulated to account for the influence of design 
changes in areas other than those of primary 
concern for a given problem. After gaining exper- 
ience with each component problem, the goal is to 
develop a completely integrated approach to 
optimize on several design considerations simul- 
taneously. Based on experience with the individual 
optimization problems, it will be possible to 
better formulate an integrated and efficient 
overall approach. Furthermore, experience will be 
gained as to the design variables having the 
largest impact on each individual problem, the 
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tradeoffs to be expected between various design 
considerations, and the capability to meet 
specified design criteria for a given problem. 

Figure 9 also shows a few potential design 
parameters that might be used for the optimisation 
problems. These include: blade geometrical prop- 
erties, primarily associated with aerodynamic 
performance; material properties, generally 
associated with blade stress; and structural prop- 
erties, associated with vibration and stress. In 
this section, only the vibration problem is consid- 
ered. The design parameters used are mass and 
bending stiffness distributions along the blade. 
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Fig. 9. Overall Approach for Rotor Blade Design 
Opt imizat ion 


Figure 10 outlines the approach to be used for 
helicopter vibration. In order to achieve the 
computer efficiency required of any useful design 
optimization tool, a simplified vibration analysis 
is used in the primary or inner loop to develop che 
vibration parameters and other criteria to be 
optimized. Since this simplified analysis may be 
performed many times faster than the forced 
response analysis, the potential savings in time is 
significant for the many iterations that may be 
required by any constrained optimization program. 
The forced response analysis G400 is then used to 
verify the vibration characteristics of the new 
blade design in the outer optimization loop, where 
closed-loop optimization can also be performed. 



Fig. 10. Optimization Approach for Vibration 


The simplified vibration analysis is based on 
the assumption that appropriate modal parameters 
can be defined that indirectly relate changes in 
vibration characteristics to changes in blade 


design. Traditionally, frequency placement has 
been used to minimize vibration response. As 
discussed above, modal shaping has also been shown 
to be important, and associated parameters have 
been identified as indicators of rotor blade vibra- 
tion characteristics. In particular, it has been 
predicted tnat vibration can be reduced by mini- 
mizing certain weighted modal integrals. In 
studies performed to date, polynomial approxima- 
tions to the airload distribution have been used as 
weightings in the modal integrals. As shown 
conceptually in Fig. 11, mode shaping is accom- 
plished by driving these generalized airloads to 
zero to desensitize the blade to vibratory air- 
loading. 



• DESENSITIZE BLADE TO VIBRATORY AIRLOADS 
BY SHAPING OF CRITICAL MODES: 


GENERALIZED _ f . riF _ 
AIRLOAD = J*N dx dx 


- 0 


Fig. 11. Modal Shaping Design Concept 

Figure 12 shows a more detailed schematic of 
the inner optimization loop presented in Fig. 10. 
A blade eigensolut ion analysis (E 159 ) is used to 
calculate blade natural frequencies and mode shapes 
for a given set of design variables. This informa- 
tion and the assumed airload distributions are used 
to calculate the appropriate modal integrals and 
the difference between the actual and optimum modal 
frequencies. Frequency placement and modal shaping 
are accomplished by simultaneously driving these 
parameters to zero via minimization of a quadratic 
performance index that consists of the weighted sum 
of the squares of each vibration parameter. The 
weighting matrix, Wz, is used to reflect the 
relative importance of each vibration parameter. 


The constrained optimization program used for 
the results presented in this paper is COPES/CONMIN 
(Refs. 4 and S), which is based on the Method of 
Feasible Directions. This program minimizes the 
performance index in an iterative manner. At each 
step, it attempts to satisfy all specified con- 
straints, which may be either explicit or implicit 
functions of the design variables. As shown by the 
dashed line in Fig. 12, blade frequencies and modal 
integrals can also be included as constraints 
rather than added to the performance index. Based 
on gradient and functional information for the 
objective and constraint functions, COPES/CONMIN 
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Fig. 12. Automated Design Procedure for Vibration 
Based on Simplified Approach 

calculates necessar*/ changes in the design vari- 
ables to further reduce the performance index at 
the current iteration. The necessary gradients are 
calculated by finite differences. 

Analytical Results 

The simplified approach outlined in Figs. 10 
through 12 was used to optimize the previously 
described baseline articulated rotor operating at a 
steady 160 kt flight condition. Thirty (30) design 
variables were used to tailor three selected modes. 
The design variables consisted of the flatwise and 
edgewise bending stiffnesses and the mass at each 
of ten (10) spanwise blade stations. Optimum 
frequency and modal integral values were specified 
for three selected modes to give a total of six (6) 
vibration parameters to be reduced. As discussed 
in the last section on blade, design for reduced 
vibration, G400 was used to identify the key modes 
and the associated frequency placement and modal 
shaping criteria for the articulated rotor 
investigated. The three modes selected were the 
first and second elastic flatwise modes and the 
first elastic edgewise mode. These modes were 
selected to reduce the response of the two inplane 
hub shears, which were identified as primary 
contributors to fuselage vibration in this 
rotorcraft. The modal integral used for each mode 
included a cubic weighting function of the blade 
spanwise loca* ! < (x) to approximate the airloading 
for this high speed flight condition. 

Analytical studies with G400 indicated the 
potential for reduced vib. 'tion response in this 
articulated rotor if the first elastic flatwise and 
edgewise frequencies could be tuned to the range of 
3.2 to 3.5/rev and 5.5 to 5.7/rev, respectively. 
These studies also showed shaping of the first and 
second elastic flatwise modes to be of prime impor- 
tance. Thus, the overall objectives of the design 
problem were to meet the specified frequency 
criteria, to drive the first two flatwise modal 
integrals with cubic weighting to zero, and to 
maintain about the same blade weight if possible. 
The second flatwise frequency and the first edge- 
wise modal integral were monitored but not included 
in the performance index. 


The automated design procedure was used to 
achieve an optimized design (modified Design A) 
having the desired dynamic characteristics shown in 
Fig. 13. As shown in this figure, the automated 
design procedure met all the criteria specified for 
the primary vibration parameters. Furthermore, 
these significant changes in dynamic characteris- 
tics were achieved while allowing only a 4 percent 
increase in blade weight. In order to obtain this 
design, blade weight was added to the performance 
index to be minimized along with the frequency and 
modal shaping criteria. In addition to practical 
upper and lower bounds placed on each design vari- 
able, constraints were applied to the first two 
flatwise modal integrals to emphasize their import- 
ance and to ensure satisfaction of a threshold 
value of -0.005. Note that this value corresponds 
to the dashed line shown in Fig. 13 and represents 
reductions of eighty and ninety percent in the 
first two flatwise modal integrals, respectively. 
As added benefits, the second flatwise frequency 
was driven away from 5/rev and the first edgewise 
modal integral was reduced by over 95 percent. 



Fig. 13. Dynamic Characteristics of the Baseline 
and Optimized Blade Designs 


Figure 14 compares the final mass distribution 
for the optimized blade design (modified Design A) 
to the distribution for the baseline production 
blade. The cross-hatched region represents blade 
root-end hardware which was not modified. While 
the blade weight for both the production blade and 
the optimized design were about the same, the mass 
distributions were significantly different. The 
automated design procedure’ shifted almost 15 lbs 
from mid-span to the outer 25 percent of the blade. 
This was required to achieve the substantial 
increase specified in the first flatwise natural 
frequency. Aj an added benefit, the increased mass 
outboard also improves rotor auto-rotation charac- 
teristics. About a 40 percent increase in edge- 
wise stiffness across most of the blade span was 
required to achieve the high frequency specified 
for the first elastic edgewise mode. Changes made 
by the analysis in flatwise structural stiffness 
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Fig. 14. Comparison of Mass Distributions for 
Baseline and Optimized Blade Designs 


along the blade span were insignificant, since 
flatwise frequency and mode shape requirements were 
accomplished through changes in the spanwise mass 
distribution. 

Tne vibration characteristics for the 
optimized blade design determined by the simplified 
automated design procedure were verified in G400. 
These characteristics were presented in Figs. 2 and 
3 and were discussed in the section on blade design 
for reduced vibration. In summary, the two inplane 
hub shear components were reduced by over 65 
percent and the vertical shear by 20 percent. As a 
result, predicted reductions in vibration on the 
order of 50 percent were achieved in cockpit and 
cabin vibration. 

Furthermore, significant reductions in flat- 
wise and edgewise bending stresses and in torsional 
stress all along the blade span were predicted. 
Reductions of nearly 50 percent were achieved at 
all the critical stress areas (outboard flatwise, 
roidspan edgewise, and inboard torsional) despite 
the lack of stress constraints and stress terms in 
the performance index. 

Closed-Loop Higher Harmonic Control 

The use of .a self-adaptive controller to 
implement higher harmonic control (HHC) in closed- 
loop fashion potentially allows significant vibra- 
tion reduction to be achieved throughout the flight 
envelope. In this approach, higher harmonic blade 
root pitch, which can be input through the standard 
swashplate configuration, is used to modify blade 
airloads and reduce harmonic blade forcing of the 


fuselage. Reference 20 presents an excellent 
review of past analytical and experimental work in 
helicopter higher harmonic control. More recently, 
the concept of closed- loop HHC has been success- 
fully demonstrated in flight tests (Ref. 21). 

In recent years. UTRC has focused on the 
analytical development, evaluation, and refinement 
of closed-loop self-adaptive higher harmonic con- 
trol algorithms. References 22 and 23 present the 
results for a numerical simulation of a closed-loop 
deterministic control algorithm. The simulation 
was based upon a Black Hawk (UH-60) aircraft flying 
at various steady flight conditions. References 7 
and 8 present the results of a more recent analyt- 
ical study involving a simulation of the H-34 rotor 
mounted on the NASA Ames Rotor Test Apparatus (RTA) 
of the 40 x 80 ft wind tunnel. This investigation 
involved the refinement and evaluation of alterna- 
tive controller configurations in order to compare 
their performance and to more fully understand the 
effects of tuning parameters within the algorithms. 
A generic controller computer code was developed to 
give the capability to readily define many dii f er- 
ent algorithms by selecting from three control 
approaches (deterministic, cautious, and dual), two 
linear system models (local and global), and sever- 
al methods of limiting control. The generic con- 
troller is currently being used in analytical stud- 
ies in preparation for open- and closed-loop flight 
tests of the Sikorsicy S-76 (Ref. 24). An overview 
of the generic controller, analytical simulation of 
closed-loop control, and results presented in de- 
tail in Ref. 7 is given below. 

Generic Active Controller 


Figure 15 shows the computer simulation used 
to evaluate and compare the performance of the 
alternative algorithms included in the generic 
active controller. This simulation of closed-loop 
control is achieved by linking the generic control- 
ler to a nonlinear aeroelastic analysis (G400), 
which simulates the rotorcraft by calculating the 



Fig, 15 Simulation of Active Vibration Control 
System 
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vibration response at a set of fixed-system sensor 
locations. Based on this response and on-line 
identification of system parameters, the active 
controller calculates and commands the HHC inputs 
required to further reduce vibration in the fuse- 
lage. These commanded inputs are imposed on the 
rotor by the G400 analysis with the inputs pre- 
scribed by the rotating system components (i.e., 3, 
4 and 5 per rev). For computat ional efficiency, a 
constant inflow model has been used in G400 . No 
measurement noise was simulated in this initial 
investigation. 

Regardless of the control approach or system 
model implemented by the generic controller, there 
are two fundamental characteristics of the active 
controller: (1) a quasi-static linear transfer 
matrix (T-matrix) relationship between the vibra- 
tion response and the HHC inputs is assumed and (2) 
the T-matrix is identified on-line to account for 
changes due to system nonlinearities or variations 
in flight condition. Each control approach can be 
based on one of two system models. The local model 
linearizes the system T-matrix about the current 
vibration response, Z, while the global modal 
linearizes about the uncontrolled vibration level, 
Z 0 (zero HHC), which must also be identified. 

As shown in Fig. 15, accurate identification 
of the T-matrix, as well as Z Q for the global 
model, is important for good vibration reduction, 
since the minimum variance control algorithms all 
depend explicitly on the estimates of these param- 
eters. The method used for estimating and 
tracking these system parameters is discussed in 
detail in Ref. 7. In short, each row of the T- 
matrix is considered to be a time-varying state 
vector, which is tracked by a Kalman filter identi- 
fication algorithm. 

Once system identification is completed, the 
required change in control for minimum vibration in 
the ith control update is calculated by a minimum 
variance control algorithm (Ref. 7). This algo- 
rithm is based on minimization of a quadratic 
performance index that consists of a weighted sum 
of the squares of the input and output variables. 
The performance index can be written in matrix 
notation as follows: 

J ’ Z i Tw Z Z i +T i (e ' X ‘ P i^ W Zjj )Y i +6 i W 0 e i +i6 i W t0 Ae i (1) 

j 

where Yj*A 0^ for the local model and ^"(sT 1)^ for 
the global model. The index 3 acts as a switching 
function dependent on the control approach used. 
The performance index J is a function of not only 
the vector of computed harmonics of vibration (Z), 
but also the vector of pitch control inputs (6) and 
incremental change in control (A0). In the first 
term, Wg is a diagonal weighting matrix used to 
reflect the relative importance of each vibration 
component. This term, referred to later as the 
vibration index, is indicative of overall effec- 
tiveness in reducing vibration. The second term is 
used to modify the controller algorithms to account 


for uncertainties in identified system parameters 
according to the underlying assumptions of the 
control approach being used. These uncertainties 
are reflected in P^,' the covariance matrix calcu- 
lated by the Kalman filter identification algo- 
rithm. The effect of this stochastic control term 
is determined by 3, and the arbitrary stochastic 
control constant X. Finally, the last two terms 
use diagonal weighting matrices Wg and W^g to 
inhibit excessive control amplitudes and rates of 
change in cont’ 1, respectively. This "internal 
limiting" is _’d nut only to satisfy hardware 
requirements, but also to enhance controller 
performance. 

For the deterministic control approach, 3 is 
set to zero, since all system parameters are 
assumed to be explicitly known despite the fact 
that only estimates for the T-matrix are available. 
In the cautious approach, it is recognized that 
some of the system parameters are only estimates, 
and control inputs are implemented more cautiously 
than for the deterministic approach. This is 

accomplished by setting 3 equal to one. The 
resulting positive stochastic control term has a 
similar effect to that of W^ 0 or W 0> depending on 
the system model, but is dependent on the uncer- 
tainty in the identified T-matrix, as reflected by 
P^. In the dual control approach, an attempt is 
made to improve long term system identification by 
actively probing the system, while maintaining good 
control. In the generic controller, this is 

achieved with a negative value for 3. The effect, 
analagous to reductions in weighting placed on 
control inputs, causes the system probing inherent 
to the dual controller. Whereas the cautious 
controller penalizes control when identification is 
poor, the dual controller increases control. 

Finally, Fig. 15 shows that the active 
controller externally limits the optimum control 
inputs calculated by the minimum variance control 
algorithm before implementing new inputs in the 
rotorcraft simulation. This is referred to as 
external limiting since it is done outside the 
minimum variance control algorithm and without 
regard to optimality of the resulting solution. 
With external limiting, satisfaction of absolute 
control limits can be ensured. In contrast, 
internal limiting, which 13 accomplished by 
appropriate tuning of the wt i ght ing mat rices , tfg 
and W^g, takes into account the desire to inhibit 
magnitudes and rates of change of control while 
calculating the optimum solution. 

Analytical Results 

The aeroelastic simulation of the rotorcraft 
in Refs. 7 and 8 was based on a fully art eted , 
four-bladed H-34 rotor mounted on the Rotor Test 
Apparat*- (RTA), which is used to represent the 
fuselage in full scale rotor tests in the NASA Ames 
40 x 80 ft wind tunnel. Vibration response infor- 
mation was calculated at six locations in the ST A. 
These components included three orthogonal direc- 
tions (vertical, longitudinal, and lateral) and 
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were 9 pread throughout the RTA (nose, main cross- 
beam, and tail). 

A steady level-flight condition was selected 
for the initial tuning and evaluation of all 
primary controller configurations. This baseline 
condition had a forward velocity of 150 kt and a 
nominal value of 0.058 for Gp/o (8250 lb thrust). 
Based on these results, a representative controller 
configuration was selected and tuned for each of 
the three control approaches. All three control- 
lers were based on the global system model. The 
deterministic controller used internal weighting on 
A0 with W^q to maintain an acceptable rate of 
change of control. The cautious controller used 
neither external nor internal A0 limiting, but 
inherently slowed down the implementation of 
control via the stochastic control term. The dual 
controller used external rate limiting to allow the 
i herent perturbations in control inputs to occur 
without excessively compromising short term con- 
trol. The performance of these controllers were 
subsequently evaluated at several steady flight 
conditions and during several short duration maneu- 
vers as discussed below. 

Steady Level Flight Conditions - Figure 16 
presents the G400 predicted results for each of the 
three controllers operating closed-loop at the 
baseline flight condition. The simulation included 
three revs of uncontrolled flight to allow initial 
numerical transients to die out before activating 
each controller at rev 4. Figure 16 jhows pre- 
dicted time histories of the vibration index J 2 and 
the amplitude of the 3 per rev HHC input commanded 
by each control approach. While not shown, 4 and 5 
per rev inputs commanded by each controller had 
similar time histories. Since the vibration index 
is a weighted sum of the squares of all the vibra- 
tion components being actively controlled, it is a 
good indicator of overall controller performance in 
reducing vibration. Note that the vibration index 
plotted involves only the first term shown in Eq. 
(1). While the other terms are important to 
overall controller performance and stability, they 
are not indicative of vibration reduction achieved 
by the controller. 

Figure 16 shows that all three controllers did 
an excellent job of reaching a new steady vibration 
level that is greatly reduced from the uncontrolled 
vibration level at rev 4. After only two revs of 
active control, both the deterministic and cautious 
controllers achieved and maintained at least a 90 
percent reduction in the vibration index. The dual 
controller required about 5 revs of active control 
to achieve the same level. By rev 10, all three 
controllers had essentially converged to a value of 
the vibration index that was only 3 percent of the 
uncontrolled value. 

This figure also shows the time history of 
3 per rev HHC amplitude commanded by the three 
controllers. The deterministic and cautious con- 
trollers smoothly increased the amplitude of all 
three control inputs, while continually reducing 




Fig. 16. Time History of Vibration Index and 3 Per 
Rev Control at Baseline Flight Condition 
(V « 150 kt, C T /0 - 0.058) 

the vibration index. In contrast, the dual con- 
troller exhibited a tendency to probe the system by 
perturbating the higher harmonic cyclic inputs. 
This probing initially resulted in a slight degra- 
dation in short term control as can be seen in the 
vibration index. After identification improved, 
system probing diminished and the final controller 
solution was as good as that of the deterministic 
and cautious controllers. The dual controller^ 
tendency to probe the system was somewhat inhibited 
by an application of external rate limits. Without 
these limits, the perturbation in control inputs 
used to probe the system were much larger and 
resulted in m \ch worse short term control before 
converging to * final solution. 

Figure 17 compares the uncontrolled 4 per rev 
vibration levels at rev 4 to those at rev 30 with 
active control. All three controllers substan- 
tially reduced vibration at all locations except 
the two lateral components that had very low 
initial levels of vibration, which were maintained. 
Reductions in vibration for the four primary compo- 
nents were between 75 and 95 percent. 

Also shown in Fig. 17 are the fixed system hub 
vibrations. Note that angular accelerations have 
been multiplied by 1 ft to be plotted in g's in 
this figure. The two largest contributors (verti- 
cal and longitudinal) were reduced oy all three 
controllers. A substantial 75 percent decrease in 
the longitudinal component was achieved, while a 
more modest 20 percent reduction was achieved in 
the vertical component. The other four components, 
which were smaller initially, remained at about the 
same levels. The reductions in vibration in the 
RTA were achieved by a combination of reduced 
forcing at the rotor hub and vectorial cancella- 
tions of hub component contributions. 
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Fig. 17. Effect of Active Control on 4 Per Rev 
Vibration at Baseline Flight Condition 
(V - 150 kt, C T /o - 0.058) 

In addition to the baseline flight condition, 
each of the three controllers was also evaluated 
over a range of forward velocities from 112 to 150 
kt at a nominal value of 0.058 for Cj/o and over 
a range of rotor thrusts having values of 0.058 to 
0.085 for Gj Jo at a velocity of 150 kt. Between 75 
to 95 percent reductions in vibration were achieved 
by all three controllers at all steady flight 
conditions. These reductions in the individual 
components in vibration correspond t-> at least a 97 
percent reduction in the vibration index at all 
steady flight conditions. Convergence to an 
acceptable control solution occurred quickly and 
smoothly. After 5 control updates, at least an 

eighty percent reduction in the vibration index was 
achieved and maintained. These results were 

obtained at all steady flight conditions with no 
retuning of the controller and with the same 
initial T-matrix developed and used at the baseline 
(V*150 kt, Cflo m Q .058) flight condition. The 
required amplitudes of 3, 4, and 5 per rev control 
increased with both velocity and rotor thrust, but 
were each less than 1.0 degree for all steady 
flight conditions. 

All the steady flight results presented above 
for the global system model are generally applic- 
able to the local model as well. It was not until 
controller performance was evaluated during the 
short duration maneuvers discussed below that any 
significant difference in controller behavior due 
to system model was noticed. From the steady 


flight conditions studied to date, no distinct 
advantage In terms of controller performance has 
been identified for either the deterministic or 
cautious control approaches. The dual controller, 
while equally effective in reducing vibration, 
tended to have worse short term control and 
somewhat more oscillatory behavior due to system 
probing. The baseline deterministic and cautious 
controllers were relatively insensitive to less 
than optimum tuning of internal parameters. 
However, it should be noted that the use of 
internal limiting in the deterministic approach 
dramatically improved controller stability and 
performance compared to that achieved with external 
limiting. The use of properly tuned internal 
weightings on control inputs significantly improved 
the deterministic controller performance according 
to all criteria: much greater vibration reduction 
in the first step of active control; faster conver- 
gence; significantly greater reduction in vibration 
at convergence; and smaller final control inputs. 
The dual controller was very sensitive to the 
tuning of \ . It remains to evaluate the effect of 
measurement noise on the performance of each 
control approach. 

Short Duration Maneuvers - Each of the three 
controllers was also evaluated during several short 
duration maneuvers, while using the same initial T- 
matrix and tuning developed at the steady baseline 
condition. Each of the maneuvers involved an 
increase in rotor thrust from the initial steady 
baseline condition, C^/o * 0.058, via step and 
ramp changes in collective pitch during an other- 
wise steady flight condition at 150 kt . These 
changes in collective pitch resulted in 40 to 50 
percent increases in rotor thrust relative to the 
baseline condition. After all transients subsided, 
the final flight condition corresponded to one of 
the steady flight conditions investigated (i.e., 
Cf/o m 0.08 or 0.085). For each of those maneu- 
vers, the active controllers remained stable, main- 
tained peak vibration response well below uncon- 
trolled levels, and reduced vibration to the same 
levels achieved at equivalent steady conditions. 
Retuning of the controllers was necessary to 
achieve satisfactory performance during some maneu- 
vers. Without retuning, the local model was much 
more oscillatory and required more time to converge 
than the global model during maneuvers. This may 
indicate that the local model is more sensitive to 
tuning at different flight conditions or perhaps 
more sensitive to inaccurate vibration response 
information due to transient effects. Detailed 
results for the various steady flight conditions 
and short duration maneuver^ are presented in Ref. 
7. The results for the maneuvers investigated 
indicate the need for further evaluation during 
extended continuous maneuvers. 

Blade Stresses and Rotor Performance - 
Increases in rotor blade stresses were noted at 
most flight conditions investigated. However, 
results al?^ suggest that the penalty of increased 
vibratory blade loads may be reduced by tailoring 
of HHC inputs with unequal Wq weighting. It may 
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also be possible to alleviate the increases in 
stress, without compromising vibration reduction, 
by including appropriately weighted terms represen- 
tative of blade stresses in the performance index 
shown in Eq. (1). While these approaches were not 
pursued, certain results did indicate that they 
might be feasible. For example, multiple control 
solutions resulting in similar vibratio.* reduc- 
tions, but having different effects on rotor blade 
stresses, were obtained. One s\ch solution 
involved an arbitrary elimination of 5 per rev 
control at the highest thrust condition investi- 
gated. The result was relatively small increases 
in blade stresses and excellent reductions in 
vibration. The increase in blade stresses were 
much smaller than those where 5 per rev control was 
implemented, even though vibration reductions for 
both cases were comparable. 

A degradation in rotor performance was also 
noted at many flight conditions. At the baseline 
condition, the application of HHC caused an 
increase in required torque on the order of 5 
percent. An analytical study of the effect of 
closed-loop HHC on rotor performance should be 
performed when using a variable iriilo;* model that 
includes unsteady aerodynamic effects. Again, it 
may be possible, if necessary, to guide the 
controller tc better control solutions in terms of 
rotor performance, as well as vibration, by 
including a term indicative of rotor torque in the 
performance index. 

Blade Appended Aeroelastic Device 

In addition to alleviating vibration through 
blade design and active pitch control, the use of 
blade appended aeroelastic devices has been analyt- 
ically explored at UTRC. One such device, as 
reported in Ref. 9, is a passive tuned tab sh jwn in 
Fig. 18. The objective of this tab is to create 
harmonic air loading of favorable amplitude and 
phase to cancel the inherent harmonic airloading 
which acts as a source of main rotor vibration. 
Physically, the passive blade tab is appended near 
the trailing edge of a standard rotor blade by some 
hinge configuration so that the tab can deflect 
freely about the hinge. The hinge could be mechan- 
ical in nature with bearings or it could be made of 
a composite material that has a large allowable 
strain such that the tab is actually "taped" to the 
blade by the composite hinge. The latitude in 
selecting the spring rate of the tab would provide 
a dynamic tuning capability; the spring rate could 
be provided either m« hanically or by the elastic- 
ity of the material for a composite hinge. 

The basis of the concept, as described in Ref. 
9, is as follows: when a rotor blade tab deflects, 
it creates an incremental airload and pitching 
moment on the rotor blade as a result of the 
increased .jmber. The pitching moment also creates 
an additional airloading on the rotor blade by 
elastic twisting to create an incremental angle-of- 
attack. The importance of this source of airload- 
ing is closely tied to the blade torsional atiff- 



Fig, 18. Blade Vibration Control Device — Passive 
Tuned Tab 


ness and natural frequency, and is secondary to 
that obtained from the effective camber change for 
this concept. When the tab deflects harmonically, 
the airloads and pitching moment created by the tab 
deflection are also harmonic. Therefore, to derive 
benefit from the tab, the tab motion must be 
correctly phased to cancel the inherent harmonic 
airloading that excites the blade flatwise modes 
and produces vibration. 

The driving forces on the tab are its *>wn 
inertial loading as the blade flaps and pitches 
(rigid body and flexible motions) and the aero- 
dynamic forces arising from blade and tab motion. 
By increasing the offset of the tab center of grav- 
ity from the hinge, the inertial forcing can be 
increased. F^r a tab located at the blade tip, 

most of the vertical harmonic motion would come 
from the response of the flexible flatwise modes 
and nearly all of the torsion motion would be due 
to the response of the blade first torsion mode. 
Hence, there is a direct relationship between the 
motion that is inertially forcing the tab to de- 
flect and the vibration that is a result of that 
same motion. Therefore, the success of this con- 
cept depends on correctly sizing and placing the 
tab along the rotor blade span and choosing its 
mass and natural frequency to achieve the maximum 
vectorial cancellation of inherent harmonic air- 
loading. 

An initial analytical evaluation of this con- 
cept has been conducted (Ref. 9) using the GAOO 
rotor aeroelastic analysis (Ref. 6), together with 
data for a realistic helicopter rotor blade 
(UH-60A, Blawk Hawk) in high speed flight (175 
kts). Variations in tab mass, frequency, and cen- 
ter-of-gravity location were investigated for two 
tab spanwise locations. While some modest reduc- 

tions in the inplane components (longitudinal and 
lateral) of vibratory hub shear were predicted, 
unacceptable increases in the vertical component 
have been predicted for the blade/tab as configured 
in the study. The reason for the vertical shear 
increase remains to be determined before any 
further investigation to determine if other config- 
urations have potential for overall vibration re- 
duction. 

Vibratory Airloads 

An ongoing activity at UTRC has been directed 
toward the development of helicopter airload meth- 


i f 




355 


odology and the advancement of the understanding of 
vibratory airloads of the rotor and fuselage. A 
representative sampling of results of this activity 
are presented in Reis. 10 to 19. This activity has 
included investigations of rotor wake and airflow, 
unsteady airloads related to dynamic stall, blade 
airload coupling with blade response and aeroelas- 
tic flexibility, and both rotor/fuselage and 
rotor/empennage interactional aerodynamics. 

Wake Induced Blade Airloads 

A recent investigation of helicopter rotor 
v/ake geometry and its influence in forward flight 
is reported in Refs. 10 and 11 and summarized in 
Ref. 12. This analytical investigation was con- 
ducted to generalize the wake geometry of a heli- 
copter rotor in forward flight and to demonstrate 
the influence of including wake deformation in the 
prediction of rotor airloads and performance. 
Predicted distortions of the tip vortex of each 
blade relative to the classical undistorted geom- 
etry were generalized for vortex age, blade azi- 
muth, advance ratio, thrust coefficient, rotor disc 
attitude, and number of blades based on a repre- 
sentative blade design. A computer module and 
charts (Ref. 11) were developed for approximating 
wake geometry and identifying wake boundaries and 
locations of blade-vortex passage. Predicted H-34 
rotor airloads for several rotor inflow/wake models 
were compared with test data for several flight 
conditions . 

An example of the tip vortex geometry in for- 
ward flight (30 kts ) , as predicted by the UTRC Weke 
Geometry Analysis (Ref. 13), is shown in Fig. 19. 
A sample isometric view of distorted tip vortices 
from the generalized wake model is shown in Fig. 
20. The characteristic wake distortion features 
observed from experimental results are prest i in 
the predicted tip vor x georat-ries . The forward 
and lateral sides of tbe wake are distorted toward 
the rotor relative to the undistorted wake model. 
This results in close blade-vortex passages which 
can introduce severe local azimuthal and spanwise 
gradients in blade airloads which are not predicted 
with uniform inflow and undistorted wake (rigid 
wake) analytical models. This is shown in Fig. 21 
where the predicted airloads (blade lift distribu- 
tions) based on the different inflow/wake models 
are presented in the form of surface contour plots. 
As indicated, inclusion of tip vortex deformation 
in the wake model results in increased higher har- 
monic content in the airload prediction. The out- 
board &4van< ing side of the rotor typically exhib- 
its the moit severe vibratory airload gradients 
with significant but lesser variations on the out- 
board retreating side. This is exemplified in the 
H-34 airload test data ahown in Fig. 22, taken from 
Ref. lu, as aco-jired at 48 kts (0.129 advance 
raLio) in flight test (Ref. 25). 

Typical of low speed transition conditions, a 
sharp "down-u,,' impulse is applied to the tip re- 
gion of the advancing blade and an "up-down" 
impulse is applied to the tip region of the 
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Fig. 19. Predicted Distorted Tip Vortex Geometry 



Fig. 20. Isometric View of Generalized Distorted 
Tip Vortices 

retreating blade (Fig. 22). For higher speed 
flight the vibratory airloads are concentrated at 
the outboard region of the advancing side where a 
characteristic "up-down" impulse occurs. This is 
exemplified in Fig. 23 where the airloads are shown 
for the H-34 rotor operating at 110 kts (0.29 ad- 
vance ratio), as measured in a wind tunnel (Ref. 
26). These airload characteristics have been shown 
by Hooper (Boeing-Vertol) in Ref. 27 to be surpris- 
ingly consistent for different helicopters with 
substantial differences in size, trim, and number 
of blades per rotor. In this reference, the com- 
bined higher harmonic components of airloading 
(harmonics 3-10) are plotted for several aircraft 
and, with some exceptions at high speed, these 
vibratory components are shown to be characteris- 
tically consistent. 

The ability to predict the vibratory airload- 
ing characteristics is depicted in Fig. 22 from 
Refs. 10 and 12, where the results of combining 
the Sikorsky Generalized Performance Analysis (air- 
loads), the UTRC Rotorcraft Wake Analysis (induced 
airflow), and the UTRC Wake Geometry Analysis (dis- 
torted wake) are presented for the H-34, 48 kt 
condition. The general ability to predict the 
induced airflow of the rotor wake in the vicinity 
of the blades is shown through comparisons with 
laser velocimeter and other test dfta in Refs. 15 
and 16. The influence of wake distortions on blade 
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airloads is describe 1 in these leferences to be 
related to the degree of proximity of the tip vor- 
tices to the rotor blades and the number of loca- 
tions of close blade-vortex passages. For steady 
level flight, the wake influence generally in- 
creases with decreasing advance ratios, decreasing 
disc attitude and increasing number of blades. The 


wake distortions and related airload influence are 
moat complex at low advance ratios, where th^ dif- 
ficulty of accurate wake geometry and airloads 
prediction is compounded by the larger number of 
close blade-vortex passages and by blade-vortex 
impingement due fo the mov' .. „ of the tip vortices 
above and then down through the rotor disk. 


O 








The prediction of tip vortex distortions, 
which ere in close proximity to the blades at 
localized regions of the rotor, results in a high 
sensitivity of predicted airloads to small differ- 
ences in tip vortex geometry and the theoretical 
blade-vortex aerodynamic interaction model used. 
This sensitivity is demonstrated in Refs. 10 and 
12. For example, the airload predictions in Fig. 
22 resulted from the use of an expanded vortex core 
model to simulate vortex bursting near close blade 
passages. Both higher and lower amplitudes of 
vibratory xirloading on the advancing and retreat- 
ing sides were 'redicted with other vortex core 
models and small variations in wake geometry 
(blade-vortex spacing). Although the need for 
further analytical refinement is evident, the de- 
gree to which the analysis, with a distorted wake 
model, was able to reproduce the measured airload 
distributions is encouraging and indicative that, 
with future emphasis on blade-vortex interaction 
modeling, wake methodolgy has the potential to 
provide a useful predictive tool for vibratory 
airloads . 

Wake geometry and blade-vortex intersection 
plots in Refs. 10 and 27 indicate the source of the 
impulsive type airloads on the advancing and 
retreating sides to be at least partially attrib- 
utable to close blade-vortex interaction. It is 
shown that the blade passes close to a nearly 
parallel tip vortex from a preceding blade in the 
first quadrant and the fourth quadrant of the 
rotor. At higher flight speeds, where the vibra- 
tory airloading is predominantly on the advancing 
side, the vibratory airloads are caused by a com- 
bination of blade-vortex interaction, blade mo- 
tions, and the characteristic blade "Mach tuck" 
phenomenon. This phenomenon is known to result 
from the aft movement of the blade aerodynamic 
center, due to compressibility, which results in a 
nose down pitching moment for positive lift and a 
corresponding blade torsional deformation. This 
produces a negative azimuthal gradient in blade 
lift. Blade sweep has been used to counte this 
effect in the advancing side negative lift region. 
In Ref. 27, it is hypothesized that fhe "up-down" 
impulse on the advancing side at high speer condi- 
tions is opposite to the "down-up" impulse at lower 
speeds due to the opposite direct ii .i of vorticity 
produced by the negative lift region at the tip at 
high speed. In high speed flight, vibratory air- 
loads on the retreating side can also be caused by 
dynamic stall as shown and predicted in Ref. 18, 

Airloads From Blade Motions 

The mechanism by which vibratory airloads 
couple with blade motions and modal response are 
described in Ref. 1. A primary source vibratory 
airloads in forward flight is the interharmonic 
coupling of airloads and blade response that 
increases with forward speed. As mentioned 
earlier, the phenomenon can be viewed as a cascade 
effect (Fig. 1), starting with the basic character- 
istic of helicopter rotors in forward flight. One 
per rev airloads are created by either blade flap- 


ping or cyclic pitch. The one per rev airioads 
cauae one oer rev response of the blade modea. The 
resulting one per rev motions create 2 and 3 per 
nv airlcada. These airloads in turn cause 2, 3, 
4, and 3 per rev response cf the blade modes and so 
on. Different blade modes respond more than others 
due to resonant type responses, and cause vi oration 
from the c<mbination of inertial, elastic, and 
direct airioading. Thus, significant vibratory 
airloads are created by blade !!x>tions and elastic 
response as well as wake induced, compressibility, 
and dynamic stall effects. 

Interactional RoLor-Fuselage Airloads 

In order to accurately predict coupled rotor- 
fuselage vibrations, it is necessary to consider 
the aerodynamic interaction of the individual com- 
ponents of the helicopter. In Ref. 17 i' is shown 
that the presence of the fuselage distorts the 
rotor airflow and wake causing two-per-rev as well 
as other harmonic airload excitations at the rotor. 
Examples from Ref. 17 of the calculated effect of 
the airframe pres, .nee on the rotor inflow veloci- 
ties and inner blade angle-of-attack distribution 
are shown in Fig. 24. Here, the influence of the 
fuselage on the rotor inflow velocities 'AV^) is 
shown as predicted using the Sikorsky fuselage 
panel method (WABAT) for a '•elatively low rotor 



FUSELAGE ANGLE OF ATTACK = -8° 
POSITIVE ^V N /V = UPWASH 
IF = 2 8 



Fig, 24. Predicted Body Induced Velocities at 
Rotor Plane and Body Effact on Rotor 
Angle of Attack Distribution 
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configuration (IF a 2.8). (The interference factor, 
IF, defined in Ref. 17, is representative of the 
ratio of the fuselage upper surface area to the 
rotor hub to fuselage dirtance.) The perturbation 
of the flow field at the rotor disc produces 
chai.ges in the blade local angle-of-at tack distri- 
butions which influence the blade response and 
airloads. Also in Fig. 24, the angle-of-at tack 
distribution at one radial station is presented as 
predicted using the UTRC Rotorcraft Wake Analysis 
coupled with the Sikorsky Generaliz d Rotor Perfor- 
mance Analysis. In addition to Iwo-per-rev, there 
is a significant third narraonic content in the 
increment rep;^senting the body induced effect. 
For low rotor configurations, this should be con- 
sidered in computations of the vibration spectrum. 

The influence of the rotor airflow on fuselage 
vibratory airloading must also be considered 
Experimental results have indicated that the rotor 
blades and wake can produce significant oscillating 
fuselage pressures at principally the fundamental 
blade passage frequency. These pressure pulses 
have been recognized as a fuselage vibration mech- 
anism, but this mechanism and its effect have not 
been fully investigated analytically. Methodology 
is currently being developed at UTRC to approximate 
the rotor wake deformation due to the fuselage and 
to calculate the unsteady fuselage pressures. 

Interactional Rotor - Empennage Airloads 

The kotorcraft Wake Analysis and laser veloc- 
imcter measurements have been used at UTRC to 
demonstrate the significant influence of the rotor 
wake on fluctuating flow velocities in the vicinity 
of tail surfaces (Ref. 16). More recently, a com- 
puter program (RTEVA) has been developed that pre- 
dicts the unsteady airloads that are imposed on the 
empennage surfaces due to its aerodynamic interac- 
tion with the main rotor wake ief. 19). A rotor 
wake program is used to determine the position and 
the strength of blade tip vortices that pass near 
the empennage surfaces. A nonlinear lifting sur- 
face analysis is utilized to predict the aerody- 
namic loads on the empennage surfaces in the 
presence of these concentrated vortices. The non- 
linear analysis was formulated to include pertinent 
effects such as suction effects cf the interacting 
vortices and the effects of time-variant shed vor- 
ticity behind the empennage surfaces. The problem 
is solved in a stepwise manner (time-domain); that 
ia, a period corresponding to one blade passage is 
divided into a large number of time intervals and 
empennage unsteady airloads are computed at each 
time step. The output of tne Analysis consists of 
chordwise and spanwise airload distributions at 
each time step. These airload distributions are 
converted into harmonic airloads that can be 
applied to excite the tail boom in a vibration 
analysis. The results of a limited correlation 
study involving the application of the computer 
program to a full-scale helicopter stabilizer indi- 
cate an encouragingly good correlation between the 
analytical vibratory airload predictions and flight 


test data. Sample results are shown in Fig. 25 
for the 6 per rev airloads of the CH-53A horizontal 
stabilizer that are induced by the wake of the six- 
bladed rotor at a high speed condition (159 kts). 



TIME (01) DEG 

Fig. 25. Unsteady Airloads Induced by Rotor Wake 
on Empennage Lifting Surface 


Concluding Remarks 

Analytical studies have demonstrated the 
potential of various approaches to alleviate heli- 
copter vibration. The approaches described herein 
should be considered to be complementary rather 
than entirely competitive. Blade design for vibra- 
tion is the moat desirable approach if vibration 
reduction can be achieved throughout the flight 
regime without penalties in performance, stresses, 
and weight. However, it may still be necessary to 
complement design optimization with other 
approaches such as active higher harmonic control 
to achieve the industry goal of 0.1 g vibration 
levels . 

A blade design study has been conducted to 
identify important contributing components and 
establish appropriate vibration criteria. The 
availability of optimization techniques makes the 
design process more tractable and allows new 
designs that satisfy specified design criteria to 
be achieved in a much more efficient fashion than 
is possible :n the traditional design process. As 
a result, it is possible to pursue more sophisti- 
cated design criteria and to achieve designs pre- 
viously unobtainable. Since it is computationally 
inefficient to base a closed-loop design optimiza- 
tion procedure on a forced response analysis, vi- 
bration criteria based on blade modal properties 
have been developed. Results from the closed-loop 
analysis are then verified in a forced response 
analysis. Results of the blade design study for a 
high speed cruise condition indicate the potential 
for developing enhanced blade designs that can 
offer significant reductions in baseline vibration 
without resorting to special vibration alleviation 
devices, radical blade geometries, or weight penal- 
ties. Efforts to define further vibration criteria 
are underway. 
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Analytical studies have also demonstrated the 
potential for active higher harmonic control to 
reduce helicopter vibration. The results of the 
generic controller study indicate the potential for 
various control approaches to reduce vibration by 
75 to 95 percent with small blade pitch amplitudes 
at a range of moderate co high speed and thrust 
conditions. Good controller performance was also 
demonstiated during short duration maneuvers. 
However, the potential for adverse effects on blade 
stresses ard rotor performance was noted at many 
flight conditions. Analytical investigations are 
planned to determine the drivers of these effects 
and to develop methods for their alleviation. 

The results for a passive tuned tab have been 
somewhat disappointing, with increases in some 
vibration components occurring at the same time 
that reductions in other components are achieved. 
Further investigation of blade appended devices, 
such as the passive tuned tab, is required. 

Computerized analyses have been used as diag- 
nostic tools to acquire a better understanding of 
the fundamentals of vibration by tracing the prin- 
cipal sources of vibration and loads and identify- 
ing the important contributing components. These 
analyses have also been used in simulating closed- 
loop higher harmonic control and in verifying the 
results of blade design changes or the effects of 
blade appended devices. Since the analyses remain 
to be validated, qualitative results regarding 
dynamic trends and characteristics of the rotor 
should be accepted more than quantitative results. 
The predicted absolute vibration and load levels 
are approximate as are the predictions of all other 
forced response helicopter analyses in use today. 
Significant progress has been made toward estab- 
lishing a blade design procedure, and preparations 
for model tests to provide vibration data are 
underway. 

Recent advancements of aerodynamic analytical 
tools are providing an understanding of the sources 
of vibratory airloads. Although further analytical 
refinment is continuing, the degree to which the 
distorted wake model is able to predict measured 
airloads is encouraging. Indications are that, 
with future emphasis on modeling blade-vortex in- 
teraction, the combination of distorted wake and 
aeroelastic response methodologies has the poten- 
tial to provide a useful predictive tool for vibra- 
tory airloads. 
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DISCUSSION 
Paper No. 22 


ANALYSIS OF POTENTIAL HELICOPTER VIBRATION REDUCTION CONCEPTS 
Anton J. Landgrebe 
and 

Mark W. Davia 


Ed Austin, U.S. Army Applied Technology Laboratory : Jack, I was wondering on your higher 

harmonic control work, first of all, what sort of a frame time were you using for updating your 
controller, and did you try investigations of various frame times to see what was necessary? 

Landgrebe : As far as frame time, you're talking about azimuth? 

Austin : Right, the computation time for the controller. 

Landgrebe : That was basically going through a one rev-type update. There was not that much 

done in [determiningj the sensitivity. It seemed to be adequate — if we had not gotten good 
results, we would have looked at other updates. 

Andy Lemnios, Kaman Aerospace : Jack, I'm absolutely delighted to have you relook at the 

dynamics end of blade design in addition to the aerodynamics end, because I think this is 
a very fertile area for some very serious vibration reductions. As you know, we feel very 
strongly about the proper inertial coupling among the various modes, and by doing so you can 
very strongly influence what happens to the hub shears and hub moments. Also we feel that 
the first flapwise frequency being above 3 per rev as we talked about last night, is a very 
Important parameter to design for. I did have a couple of questions for you. In your look 
at loads being transmitted to the fuselage from the rotor, did you also include the pitch link 
loads in that analysis? 

Landgrebe : No, that was not included. 

Lemnios : Typically, I know that in some helicopters, for example, those pitch link loads can 

account for 25 to 30J of the excitation forces. They're very significant and they cannot be 
ignored. 

Landgrebe : Yes, we recognize that. They have been looked at in other studies as far as our 

unsteady aerodynamics studies, and you're right, they can be quite formidable. 

Lemnios : The second question I had was, in this particular model did you have a lead-lag damper 

on this or not, or was it a bearingless rotor that you were looking at? Do you remember? 

Landgrebe : In which study, in the vibrations study? 

Lemnios : The /it rat ions s l ;udy, yes. 

Landgrebe : That was basically an S-76 type articulated blade. 

Lemnios : The reason I ask is, did you put in the appropriate frequency characteristics of 

the lag damper? I know from our experience, lag dampers typically have a very uniform damping 
characteristic, and essentially zero spring rate, at the low frequencies; but that lag 
characteristic falls off, and the spring rate builds up to a liquid spring at the 4 aiw »r 

rev frequency and by ignoring the 5 per rev frequency, again you may be overlooking som y 

significant vibratory forces. 

Landgrebe : Bob [Taylor] can probably answer that question better than I can, but the point 

is that in the G400 analysis, there is a representation for the lag damper, but it doesn't take 
into account all the features that you have mentioned. 

Bob Taylor, Boeing Vertol : I personally haven't seen any Sikorsky data or any other data at 

high speed where the 5 per rev vibration responses are important. 

Lemnios : Strictly from our own experience, again, our frequency inplane Is on the order of 

I think it was 5.4 per rev with the lag damper characteristics thrown in there, and that could 
be a significant contributor. 

Taylor : I do believe that the [G400] model Just used a linear lag damper. 

Lemnios : Last but not least, the conclusion on your passive tuned damper may be appropriate, 

but you ought to also think about our favorite, which is a controlled tab instead of a passive 
tuned tab. 
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Peretz Filedmann, University of California, Los Angeles : I like the way you do the optimiza- 

tion with that modal -shaping type of concept. One of the things I was wonder ing--when you 
change the mass and stiffness distribution of your blade in order to reduce your vibration 
levels, you probably affect also the aeroeiastic stability of the blades. Yet among the various 
constraints which you have enforced and listed in your slide, you had no aeroeiastic stability 
constraint. 

Landgrebe : That i3 correct, Peretz, that was not included within the constraint itself. 

The tack was t^ken to run the analysis and look to see if the analysis indicated any stability 
problems, which it did not. So that was not within the [constraint], but it could very well be 
put in in the future. 

Friedmann : You miss the point of what the constraint means. It means that if you have a 

certain given aeroeiastic stability margin, as a result of going through your procedure, that 
margin is not diminished. You still have a stable blade, yet you wouldn’t like to trade off 
aeroeiastic stability for vibration reduction of if you want to do it you should put some 
penalty on it. 

Landgrebe : Obviously the concentration here was on vibration. I showed in one of the slides 

that really what you want to do is include performance, stability, loads, put them all in at the 
same time and that is the direction to go in, but we are not to that point yet. 
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